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ABSTRACT 

m 
o 

^ . We consider the evolution of a spherically expanding plasma cloud, where 

5-1 , there is continuous injection of non-thermal electrons. We compute the time 

^ . dependent electron distribution and resultant photon spectra taking into account 

^ ! synchrotron, adiabatic and inverse Compton cooling. This model is different 
from previous works where, instead of a continuous injection of particles, a short 

. injection period was assumed. 

! We apply this model to the radio/optical knots in the large scale jets of AGN, 

! detected in X-rays by Chandra and find that the overall broadband spectral fea- 

^ ! tures can be reproduced. It is shown that for some sources, constraints on the 

'. X-ray spectral index (by a longer Chandra observation) will be able to differen- 

O 



6 



tiate between the different models. This in turn will put a strong constraint on 
the acceleration mechanism active in these sources. 

Subject headings: Galaxies: active - galaxies: jets - X-rays: galaxies 



^ . 1. Introduction 



The detection of Kilo-parsec (kpc) scale jets with knots in several active galactic nuclei 
(AGN) by the C/iandra observatory (Chartas et al. 2000; Tavecchio et al. 2000; Sambruna et al. 2002) 
has opened a new window on the nature of these phenomena. The X-ray emission from 
jets could be due to Comptonization of the optical/radio synchrotron photons i.e., syn- 
chrotron self-Compton (SSC), Comptonization of cosmic micro- wave background photons 
(IC/CMB), or an extension of the synchrotron radio/optical emission. Earlier observa- 
tions by the Einstein and ROSAT satellites, having comparatively limited resolution and 
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sensitivity, had not been generally able to differentiate amongst these mechanism (e.g. 
(Biretta et al. 1991; Harris & Stern 1987)). Based on the recent Chandra observation, it 
has been argued that the SSC interpretation would require large jet powers and magnetic 
fields much lower than the cquipartition values (Tavecchio et al. 2000; Schwartz et al. 2000). 
A plausible attractive alternative is to interpret the X-ray emission as being due to IC/CMB 
or direct synchrotron, in which case a significantly smaller jet power and near-equipartition 
magnetic fields are required (Celotti et al. 2001; Pesce et al. 2001; Sambruna et al. 2002). 
For knots where the X-ray flux is greater than the extrapolation of the radio/optical spec- 
tra to X-ray wavelengths, the emission cannot be due to synchrotron emission and hence, 
the IC/CMB model is favored. On the other hand, for those knots whose X-ray flux is 
lower than this extrapolation, it is also possible to interpret the emission as being due to 
synchrotron from a non-thermal electron distribution with a high energy cut off. In fact, 
it has been argued (Pesce et al. 2001) that for the knots of 3C271, this is a more plau- 
sible explanation, since the IC/CMB model requires for this source an exceptionally large 
Doppler factor. In low-power FRI jets also, synchrotron origin of X-rays is currently accepted 
(Worrall et al. 2001). However, for the gigahertz-peaked spectrum (GPS) radio source PKS 
1127-145, which is suspected to be a young FRI radio galaxy, the X-ray emission is probably 
due to IC/CMB (Siemiginowska et al. 2002). 

The interpretation of X-ray emission for some sources as being due to synchrotron 
emission is interesting, since in this case the X-rays are produced by particles having much 
higher energies (by a factor of «i 10^) than those emitting in the radio. The presence of 
high energy electrons puts constraints on the age of the knot: since these electrons cool 
efficiently, they are depleted in time unless they are replenished. The depletion leads to a 
time-dependent high energy cutoff in the non-thermal electron distribution which reflects as 
a high frequency exponential cutoff in the observed spectra. Since the observed -ray flux 
is less than the extrapolation of radio/optical flux in these sources, this exponential high 
frequency cutoff is constrained to be in the soft X-ray regime. This can be translated to 
a high energy cutoff in the electron distribution which in turn gives an estimate of the age 
of the knot. In this model, a knot is formed when a short duration acceleration process 
enhances the non-thermal electron density in a jet. These non-thermal electrons move with 
a bulk speed v c along the jet. Thus, from the age of the knot, one can determine the 
location in the jet of the short acceleration process. The distance from the central object 
and the short duration ( << the age of the knot) naturally puts strong constraints on any 
models of the acceleration process. 

The model outlined above, may be confirmed (or ruled out) by future long duration 
Chandra observations. The model predicts that the X-ray spectra of the knots should be 
exponential (i.e. steep). The photon spectral slope measured during preliminary short 
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duration observations of 3C 371 (Pesce et al. 2001), is around F = 1.7 ± 0.4, in apparent 
contradiction to this prediction. However, longer observations are required to confirm this 
result. Moreover, the model requires the coincidence that the age of the knot be equal to 
the time required for X-ray emitting electrons to cool. A larger survey of X-ray jets have 
to be sampled to confirm whether this is statistically plausible. In this model, it is assumed 
that the acceleration time-scale is short (i.e. tacc « 10^° sec) and the possible expansion of 
the plasma is not taken into account. The acceleration mechanism of these knots is largely 
unknown, but if the acceleration is due to internal shocks (e.g. (Spada et al. 2001)), then 
the acceleration time scale tacc ~ Rs/c-i where Rg is a typical shell size. If Rs ~ kpc, then 
tacc could be as large as 10^^ sees. 

In this work, we present an alternative model where the knot is assumed to be a uniform 
expanding sphere with continuous injection of non-thermal particles. The time-dependent 
electron distribution and the resultant spectra are computed taking into account synchrotron, 
IC/CMB and adiabatic cooling due to the expansion of the sphere. Since there is continuous 
injection of particles, there is no high energy cut-off but, instead, there is a time dependent 
break in the electron distribution where synchrotron/adiabatic cooling is important. Our 
motivation here is to show that this model also explains the observed spectra and can be 
distinguished from one time injection models, by future observational constrains on the X-ray 
spectral index. 

In the next section, the model is described and the predicted spectral energy distri- 
butions arc compared with observations. §3 summarizes and discusses the main results. 
Throughout this work. Ho = 75 km s^^ Mpc^^ and qo = 0.5 are adopted. 

2. Spherically expanding plasma model 

We consider a plasma cloud moving relativistically along the jet with a bulk Lorcntz 
factor r. In the rest frame, it is assumed that the plasma uniformly occupies an expanding 
sphere of radius R{t) = Ro + PexpCt, where i?„ is the initial size of the sphere and /3exp << 1- 
At t = 0, there are no non-thermal particles in the system. A continuous and constant 
particle injection rate [ i.e. number of particles injected per unit time] for i > is assumed, 
with a power-law distribution of energy. 



Q{'y)dj — Kj ^dj for 7 > 7, 



(1) 
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where 7 is the Lorentz factor of the electrons. The kinetic equation describing the evolution 
of the total number of non-thermal particles in the system, A'"(7, t), is 

^^^ + ^[P(7,t)iV(7,i)] = g(7) (2) 
P(7, t) is the cooling rate given by 

P{i, t) = -{^s{t) + iic{t) + iA{t)) (3) 

where 7s'(t), 7/c'(t) and 7a (^)! the cooling rates due to synchrotron, inverse Compton of the 
CMB and adiabatic cooling, respectively, are given by 

ls{t) = TT^^7' 4 
3 rUeC Stt 



i.cit) = f ^rVTl,(.)7= (5) 

T.(t) - ^ (6) 

Here, the evolving magnetic field is parameterized to be B{t) = BQ(R(t) / Rq)"^ and Tcmbi^) — 
2.73(1 + z) is the temperature of the CMB at the redshift z of the source. Note that the 
time t and other quantities in the above equations are in the rest frame of the plasma. 

In this work, Eqn (2) has been solved numerically for N{'y,t) using the technique given 
by (Chang & Cooper 1970) and the resultant synchrotron and inverse Compton spectra are 
computed at an observing time t — to- As shown by (Dermer 1995), since the CMB radiation 
is not isotropic in the rest frame of the plasma, the inverse Compton spectrum is also beamed. 
Finally the flux at the earth is computed taking into account the Doppler boosting (e.g. 
(Begelman &; Blandford 1984)), characterized by the Doppler factor S = [r(l — f3cos9)]~^, 
where (3c is the bulk velocity and 9 is the angle between the jet and the line of sight. 

While the total non-thermal particle distribution has to be computed numerically, a 
qualitative description is possible by comparing cooling timescales with the observation time 
to- The cooling time-scale due to synchrotron and inverse Compton cooling at a given time 
t and Lorentz factor 7, is tc{t,j) 7/(75 + 7/C7). Then, 7c, defined as the 7 for which this 
cooling time scale is equal to the observation time, tc(^o)7c) ~ to, becomes. 
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The adiabatic cooling time-scale also turns out to be ^ to, since ta ~ R{t)/PexpC ~ to for 
Ro << R(to). Thus, the non-thermal particle distribution at time t — to can be divided into 
three distinct regions: 

(i) 7 << 7c: In this regime, radiative coohng is not important and N{'y,to) ~ K^~Ho. The 
corresponding spectral index for both synchrotron and inverse Compton is a = (p — l)/2. 

(ii) 7 >> 7c: In this regime, either synchrotron or inverse Compton cooling is dominant 
and N{'y,to) oc 7~(?'+^). The corresponding spectral index for both synchrotron and inverse 
Compton is q; = p/2. 

(iii) 7 7c In this regime, either synchrotron or inverse Compton cooling as well as adiabatic 
cooling are important, and the spectral slope is in the range {p — l)/2 > a > p/2. 

The computed spectra depend on the following ten parameters: the observation time 
to, the magnetic field at the time of observation Bf = B{t = to), the magnetic field variation 
index m, the radius of the knot at the time of observation, Rf = R{t = to), the index p, the 
minimum Lorentz factor 7^i„, the Doppler factor S, the bulk Lorentz factor F, expansion 
velocity Pe^p ^-nd the normalization of the injection rate K. On the other hand, there are 
only three observational points namely, the radio/optical and X-ray fluxes. Clearly the 
parameters are under-constrained and it is not possible to extract meaningful quantitative 
estimates. However, the motivation here is to show that this model can explain the observed 
data with reasonable values of the above parameters. 

In Figure 1, the computed spectra are compared with the data for different knots for 
four sources. The values of the parameters used are tabulated in Table 1. The injected 
power in non-thermal particles in the rest frame is 



where Up, Ue and Ub are the energy densities of the protons, electrons and the magnetic 
field, respectively. Here, it has been assumed that the protons are cold and the number of 
protons is equal to the number of electrons. The jet power ranges from 10^^ to 2 x 10^^ ergs 
sec~^, while the injected power is generally three orders of magnitude lower. This means 
that the non-thermal acceleration process is inefficient and most of the jet power is expected 
to be carried to the lobes. The magnetic field Bf is nearly equal to the equipartition values. 

Like the results obtained by (Sambruna et al. 2002) and (Pesce et al. 2001), the X-ray 
emission from the knots in 3C371 and Knot A of 1136-135 are identified as being due to 





while the power in the jet can be approximated to be (e.g. (Celotti et al. 1997)), 

Pjet^nRT^PciUp + U. + UB) 



(9) 
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synchrotron emission. However, in this case the predicted X-ray spectral index is ax — 
aR + 1/2 instead of being exponential. Note that this relation between the spectral indicies 
is independent of the parameters used to fit the data. For the rest of the sources, like 
the earlier results (Sambruna et al. 2002), the X-ray emission is attributed to IC/CMB. 
However, for some of the sources the optical spectral index is now ao = Q^r + 1/2, instead 
of being exponential. 

3. Summary and Discussion 

In this work, it is shown that the observed radio, optical and X-ray fluxes of knots in kpc 
scale jets in AGN, can be explained within a framework of a model where there is continuous 
injection of non-thermal particles into an expanding spherical plasma. This interpretation 
can be confirmed (or ruled out) vis-a-vis one-time injection models, by future measurements 
of the radio optical ao and X-ray ax spectral indicies. In particular, the following cases 
are possible: 

(i) aR ~ ax'- In this case, the X-ray emission is probably due to IC/CMB. Both the 
continuous injection and one-time injection models are equally viable. 

(ii) ax ~ tt/j + 1/2 : In this case, the X-ray emission would be due to synchrotron emission 
from electrons in the cooling dominated region. The continuous injection scenario will be 
favored. 

(iii) ax > aR + 1/2 : In this case, when the X-ray emission is exponentially decreasing, it 
should be attributed to the high energy cutoff in the electron distribution. The one-time 
injection scenario will be favored. 

(iv) ax < oiR '■ In this case, when the X-ray emission is exponentially increasing, it should be 
attributed to the low energy cutoff (7mm) in the electron distribution and the X-ray emission 
should be due to IC/CMB. Both the continuous injection and one-time injection models are 
equally viable. 

Similar arguments can be put forth for the optical spectral index ao as compared to 
the radio. It should be noted that for some older systems the one time-injection would be 
the natural scenario, while for younger systems the continuous one would be more probable. 
The technique described above will be able to differentiate between the two and a generic 
constraint on the acceleration time-scales and typical age of the knots may be obtained. A 
generic model where the injection rate decays in time, may then be used to fit the obser- 
vations. The measurement of spectral indicies at different wave-lengths will also reduce the 
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number of unconstrained parameters in the model fitting, leading to reliable estimates of the 
system parameters. 

Constraints on the acceleration time-scale and the age of the knots would be the first 
step toward understanding the driving mechanism active in these sources. It will then be 
possible to compare these time-scales with theoretical results from analytical or numerical 
computations, which may finally lead to an understanding of the origin and nature of jets 
in AGN. 
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Table 1. Parameters for model fitting 



Source/ Knot 


B 


/ 


'Jmin 


p 


to 


S 


r 




P 


Pjet 




1136-135 A 


0, 


.9 


2.0 


2.4 


0.2 


5 


5 


0, 


,8 


44. 


,7 


47.5 


0.4 


B 


4, 


.0 


20.0 


2.9 


9 


5 


5 


0. 


,1 


44. 


,2 


47.9 


0.5 


1150+497 A 


2, 


.5 


30.0 


2.85 


9 


5 


3.5 


0. 


,1 


44. 


,2 


47.3 


0.4 


B 


4, 


.3 


30.0 


3.3 


9 


5 


3.5 


0. 


,1 


44. 


,1 


47.3 


0.75 


1354+195 A 


1, 


.7 


40.0 


3.0 


9 


3.5 


2 


0. 


,1 


45. 


,8 


48.4 


0.04 


B 


8, 


.0 


25.0 


3.2 


9 


3.5 


2 


0. 


,1 


44. 


,7 


47.5 


0.63 


3C 371 A 


1, 


.3 


10.0 


2.4 


12 


3.5 


3.5 


0. 


,1 


42. 


,8 


46.4 


0.8 


B 


1, 


.0 


10.0 


2.4 


1 


3.5 


3.5 


0. 


,5 


43. 


5 


46.0 


0.9 



Note. — Columns:- 1: source and knot name taken from (Pesce et al. 2001) 
for 3C371 and the rest from (Sambruna et al. 2002). 2: magnetic field in units 
of 10~^ Gauss at the observation time, Bf — B{t — to)- 3: minimum Lorentz 
factor 7mm- 4: non-thermal injection index p. 5: observation time to in units 
of 10^^ sec. 6;: Doppler factor S. 7: bulk Lorentz factor F. 8: log of the 
injected power in ergs s~^ Pinj- 9: log of the total jet power in ergs s~^ Pjet- 
10: ratio of the magnetic field to the equipartition value. For all cases, the 
magnetic field variation index m and the size of the source at t = to is fixed 
at 1.5 and 5 x 10^^ cm, respectively. 
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Fig. 1. — The observed fluxes in radio, optical and X-ray compared with model spectra using 
parameters given in Table 1. The data for 3C 371 are taken from (Pesce et al. 2001), while 
the rest are taken from (Sambruna et al. 2002). Triangles (dots) correspond to knot A (B). 
Errors are typically 30% or larger. 



